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£N| , Abstract 



q ■ The vacuum expectation values of the held squared and the energy-momentum tensor are 

investigated for a scalar held with Dirichlet boundary conditions and for the electromagnetic 
held inside a wedge with a coaxial cylindrical boundary. In the case of the electromagnetic 
held perfectly conducting boundary conditions are assumed on the bounding surfaces. By 
using the Abel-Plana-type formula for the series over the zeros of the Bessel function, the 
vacuum expectation values are presented in the form of the sum of two terms. The hrst one 
corresponds to the geometry without a cylindrical boundary and the second one is induced 
by the presence of the cylindrical shell. The additional vacuum forces acting on the wedge 
sides due the presence of the cylindrical boundary are evaluated and it is shown that these 
C"| | forces are attractive for both scalar and electromagnetic helds. 

PACS numbers: ll.10.Kk, 03.70. +k 

o 

^ ! 1 Introduction 



The Casimir effect has important implications on all scales, from cosmological to subnuclear, 
and has become in recent decades an increasingly popular topic in quantum field theory. Since 
the original work by Casimir [lj many theoretical and experimental works have been done on 
this problem (see, e.g., [3] and references therein). In particular, a great deal of attention 
received the investigations of quantum effects for cylindrical boundaries. In addition to tra- 
ditional problems of quantum electrodynamics under the presence of material boundaries, the 
Casimir effect for cylindrical geometries can also be important to the flux tube models of con- 
finement [3] and for determining the structure of the vacuum state in interacting field theories 
[5]. The calculation of the vacuum energy of electromagnetic field with boundary conditions 
defined on a cylinder turned out to be technically a more involved problem than the analogous 
one for a sphere. First the Casimir energy of an infinite perfectly conducting cylindrical shell 
has been calculated in Ref. [6] by introducing ultraviolet cutoff and later the corresponding 
result was derived by zeta function technique [7J (for a recent discussion of the Casimir energy 
and self-stresses in the more general case of a dielectric-diamagnetic cylinder see [8] and refer- 
ences therein). The local characteristics of the corresponding electromagnetic vacuum such as 
energy density and vacuum stresses are considered in [9] for the interior and exterior regions 
of a conducting cylindrical shell, and in [10] for the region between two coaxial shells (see also 
|llj). The vacuum forces acting on the boundaries in the geometry of two cylinders are also 
considered in Refs. [12]. The scalar Casimir densities for a single and two coaxial cylindrical 
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shells with Robin boundary conditions are investigated in Refs. [1'6\ I14|. Less symmetric con- 
figuration of two eccentric perfectly conducting cylinders is considered in [12J . Vacuum energy 
for a perfectly conducting cylinder of elliptical section is evaluated in Ref. fl5j by the mode 
summation method, using the ellipticity as a perturbation parameter. The Casimir forces acting 
on two parallel plates inside a conducting cylindrical shell are investigated in Ref. [16]. The 
Casimir effect in more complicated geometries with cylindrical boundaries is considered in |17j . 

Aside from their own theoretical and experimental interest, the problems with this type of 
boundaries are useful for testing the validity of various approximations used to deal with more 
complicated geometries. From this point of view the wedge with a coaxial cylindrical boundary is 
an interesting system, since the geometry is nontrivial and it includes two dynamical parameters, 
radius of the cylindrical shell and opening angle of the wedge. This geometry is also interesting 
from the point of view of general analysis for surface divergences in the expectation values of local 
physical observables for boundaries with discontinuities. The nonsmoothness of the boundary 
generates additional contributions to the heat kernel coefficients (see, for instance, the discussion 
in [18] and references therein). In the present paper we review the results of the investigations 
for the vacuum expectation values of the field squared and the energy-momentum tensor for the 
scalar and electromagnetic fields in the geometry of a wedge with a coaxial cylindrical boundary. 
In addition to describing the physical structure of the quantum field at a given point, the energy- 
momentum tensor acts as the source of gravity in the Einstein equations. It therefore plays an 
important role in modelling a self-consistent dynamics involving the gravitational field. Some 
most relevant investigations to the present paper are contained in Refs. [2] [T9 | [20 ]. [2T ] [22 ] [23 ] [2^] . 
where the geometry of a wedge without a cylindrical boundary is considered for a conformally 
coupled scalar and electromagnetic fields in a four dimensional spacetime. The Casimir effect 
in open geometries with edges is investigated in [25]. The total Casimir energy of a semi- 
circular infinite cylindrical shell with perfectly conducting walls is considered in [26J by using 
the zeta function technique. The Casimir energy for the wedge-arc geometry in two dimensions is 
discussed in [27] . For a scalar field with an arbitrary curvature coupling parameter the Wightman 
function, the vacuum expectation values of the field squared and the energy-momentum tensor 
in the geometry of a wedge with an arbitrary opening angle and with a cylindrical boundary 
are evaluated in [28] [29]. Note that, unlike the case of conformally coupled fields, for a general 
coupling the vacuum energy-momentum tensor is angle-dependent and diverges on the wedge 
sides. The corresponding problem for the electromagnetic field, assuming that all boundaries 
are perfectly conducting, is investigated in [30]. The scalar Casimir densities in the geometry 
of a wedge with two cylindrical boundaries are discussed in [31] . The closely related problem of 
the vacuum densities induced by a cylindrical boundary in the geometry of a cosmic string is 
investigated in Refs. [32] for scalar, electromagnetic and fermionic fields. 

We have organized the paper as follows. The next section is devoted to the evaluation 
of the Wightman function for a scalar field with a general curvature coupling inside a wedge 
with a cylindrical boundary. By using the formula for the Wightman function, in section [3] we 
evaluate the vacuum expectation values of the field squared and the energy-momentum tensor 
inside a wedge without a cylindrical boundary. The vacuum densities for a wedge with the 
cylindrical shell are considered in section HI Formulae for the shell contributions are derived and 
the corresponding surface divergences are investigated. The vacuum expectation values of the 
electric and magnetic field squared inside a wedge with a cylindrical boundary are investigated in 
section [5] assuming that all boundaries are perfectly conducting. The corresponding expectation 
values for the electromagnetic energy- momentum tensor are considered in section [6] The results 
are summarized in section [7] 



2 



2 Wightman function for a scalar field 



Consider a free scalar field <p(x) inside a wedge with the opening angle 4>o and with a cylindrical 
boundary of radius a (see figure [Q). We will use cylindrical coordinates (x 1 , x 2 , . . . , x D ) = 
(r, (f>, zi, ■ ■ ■ , zn), N = D — 2, where D is the number of spatial dimensions. The field equation 
has the form 

(V*V< + m 2 + £R) <p(x) = 0, (1) 

where R is the scalar curvature for the background spacetime and £ is the curvature coupling 
parameter. The special cases £ = and £ = £ c = (D — 1)/AD correspond to minimally and 
conformally coupled scalars respectively. 

In this section we evaluate the positive frequency Wightman function (0\ip(x)(p(x')\0) assum- 
ing that the field obeys Dirichlet boundary condition on the bounding surfaces: 



<P\ 



0=0 



4>=(f>o 



V>\i 



0. 



(2) 



The vacuum expectation value (VEV) of the energy-momentum tensor is expressed in terms of 
the Wightman function as 



(0|T ifc (x)|0> = lim V l V k (0\i P (x)i P (x')\0) + 



tow* - £ViV A 



(0|^ 2 (x)|0). (3) 



In addition, the response of a particle detector in an arbitrary state of motion is determined by 
this function. In (|3|) we have assumed that the background spacetime is flat and the term with 
the Ricci tensor is omitted. The Wightman function is presented as the mode sum 



(0|99(xMx')|0) = ^v9 Q (x)^(x'), 



(4) 



where {ip a (x), ip^(x)} is a complete orthonormal set of solutions to the field equation, satisfying 
the boundary conditions, a is a set of the corresponding quantum numbers. 




Figure 1: Geometry of a wedge with the opening angle </>o and cylindrical boundary of radius a. 



2.1 Interior region 

In the region ^ r ^ a (region I in figure [T]), the eigenfunctions satisfying the boundary 
conditions ([2]) on the wedge sides (j> = 0, 0q have the form 



<Pa(%) = fiaJqniir) sin(gn^) exp (ikrn — iut) , 
u = vV + k m = |k| 2 + m 2 , q = -n/4>Q, 



(5) 
(6) 
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where a = (n, 7, k), —00 < kj < 00, n = 1, 2, • • • , k = (k\, . . . , kjy), ry = (2^, . . . , zpj), and 
Jl(z) is the Bessel function. The normalization coefficient B a is determined from the standard 
Klein-Gordon scalar product with the integration over the region inside the wedge and is equal 
to 

The eigenvalues for the quantum number 7 are quantized by the boundary condition ([2]) on the 
cylindrical surface r = a. From this condition it follows that 

7 = X n ,j/a, j = 1,2, ••• , (8) 

where A nj - are the positive zeros of the Bessel function, Jq n (X n j) = 0, arranged in ascending 
order, X n j < X n j+i. 

Substituting the eigenfunctions © into mode sum formula (j3|) with the set of quantum 
numbers a = (n,j, k), for the positive frequency Wightman function one finds 

/oo 00 
Ae* Ar ll sm(W>) sin(gn0') £ J gn ( 7 r) J (?n ( 7 r')e-^ A *, (9) 
n=l j=l 

where 7 = X n j/a, and Ary = ry — rji, At = t — if . In order to obtain an alternative form 
for the Wightman function we apply to the sum over j a variant of the generalized Abel-Plana 
summation formula 1 33 1 



l .7/2 { . \ = / f(z)dz + -Res z=0 



^ Xn,j J^i{X n ,j) Jo 



f(z) Yqn ^ 



+- / & / , / [e-^/^e"/ 2 ) - e <?n "/(ze-" /2 )], (10) 

where Yz(z) is the Neumann function, and Ii(z), Ki(z) are the modified Bessel functions. The 
corresponding conditions for the formula (llOjl to be valid are satisfied if r + r' + \t — i!\ < 2a. 
In particular, this is the case in the coincidence limit t = t' for the region under consideration, 
r,r' < a. Formula (|10p allows to present the Wightman function in the form 



where 



(0\<p(x)<p(x')\0) = (0 w \cp(x)<p(x')\0 w ) + (ip(x)(p(x')) cyh (11) 



(0 w \cp(x)<p(x')\0 w } = - / — ^e* kAr H / dz 



00 

^ sin(qn^) sin(gn0') J ?n (zr)J qn (zr r ), (12) 



00 

x 

re=l 



and 



■ , , , / - I ^ ikArn P > ^COsh(At^ 2 -^ 
(^(X)^(x'))cyl = / 7^V e " / dZ 



7T 

00 



i^) N Jk m ' ^7^4 

^ sin(gn0) sin(gn0 / ) I qn (zr) I qn {zr') K i n ^ za ) m (13) 



n=l 
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In the limit a — > oo for fixed r,r', the term (ip(x)(p(x')) cy i vanishes whereas the part (|12|) does 
not depend on a. Hence, the term (0 w \ip(x)ip(x')\0 w ) is the Wightman function for the wedge 
without a cylindrical boundary with the corresponding vacuum state |0^). Consequently, the 
term (ip(x)ip(x')) cy \ is induced by the presence of the cylindrical boundary. For points away the 
cylindrical surface this part is finite in the coincidence limit and the renormalization is needed 
only for the part coming from the term (|12p . 



2.2 Exterior region 

In the region outside the cylindrical shell (region II in figure Q]): r > a, ^ (fr ^ (fro, the 
eigenfunctions satisfying boundary conditions ([2]) are obtained from ([5|) by the replacement 

Jqniir) -> 9qn(.7r,ja) = J qn h r ) Y qn(ia) - J qn {ia)Y gn (-fr). (14) 

Now the spectrum for the quantum number 7 is continuous and 

Pi = {2 " rDl (15) 

(j) UJ [</f n (7«) + Ygni^a)] 

Substituting the corresponding eigenfunctions into the mode sum formula (HD, the positive 
frequency Whightman function in the exterior region is presented in the form 



1 f d N k 00 
(0\<p(x)cp(x')\0) = — J j^—yj e ikAr H sm(qn(fr) sm(qn<f>') 



OCj 



d jgqn(jr,-fa)g qn (-fr',ja) e » A V7 2 +fcj 

7 J2 n ( 7 «)+n 2 n(7«) Vl 2 + kl ' 

To find the part in the Wightman function induced by the presence of the cylindrical shell, we 
subtract from (|16|) the corresponding function for the wedge without a cylindrical shell, given 
by (|12p . This allows to present the Wightman function in the form (jlip with the cylindrical 
shell induced part 



<*M*')>- = -^Jnw e "L dz 



d Nk „ikAr„ [°° ^zcoshjAWz 2 -k 2 ) 

Vz 2 - k 2 

Iqn(za) 

' K qn {za) ' 



00 

x ^2 sin(qn(fr) sm(qn(fr')K qn (zr)K qn (zr') ■ (17) 

n=l 



As we see, the expressions for the Wightman functions in the interior and exterior regions are 
related by the interchange I qn +± K qn of the modified Bessel functions. 



3 VEVs inside a wedge without a cylindrical boundary 

In this section we consider the geometry of a wedge without a cylindrical boundary. For integer 
values of q, after the explicit summation over n, the Wightman function is presented in the form 



{0 w \<p(x)<p(x')\0 v 



m 



(D-l)/2 



(2vr)( D + 1 )/ 2 



i^£(-i) i+1 £ 



q ^K {D _ 1)/2 {mJu { p 2 + |Ar 



(At)= 



3=1 



1=0 



u 



+ 



(At)* 



(£>-l)/4 



(18) 



where u\ = {r 2 + r' 2 — 2rr' cos[2ttI / q + (fr + (— 1) J (f)']} 1 / 2 . Note that the Wightman function in 
the Minkowski spacetime coincides with the term j = 1, I = in formula (1181) . 
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Taking the coincidence limit x' — ► x, for the difference of the VEVs of the field squared, 

(^ 2 )S = (0 w \<p 2 (x)\0 w ) - (0mW 2 (x)\0 m ), (19) 

where |0jv/) is the amplitude for the vacuum state in the Minkowski spacetime without bound- 
aries, we find 



m 



D-l 



/-2u»-,_ m- - A . +1 fi -g (g - 1)/2 (2mrBintf , 

{ * )ien (2*)(»+w2^ ) ^ (2mrsin ^) )( D-i) /2 • 

In this formula, the prime means that the term j = 1, I = has to be omitted, and we use the 
notation 

^) =7 rZ/g+(l + (-l)i)0/2. (21) 
For a massless field, from (|20p we find 

D-l\ 2 9-1. ,_ 

(22) 



(47r)^r«-l^^ sin^" 1 ^' 

Note that the terms in this formula with j = 2, I = and j = 2, I = q — 1 are the corresponding 
VEVs for the geometry of a single plate located at <p = and <j) = 00 5 respectively. In the case 
D = 3 for the renormalized VEV of the field square one finds [20] 

W )rJ = • ( 23 ) 

Near the wedge boundaries <j> = <p m , m = 0, 1 (<pi = 0) the main contribution in (j22H comes 
from the terms j = 2, I = and Z = qr — 1 for m = and m = 1 respectively, and the renormalized 
VEV of the field squared diverges with the leading behaviour (^jS oc \<p — (f) m \ l ~ D ■ The surface 
divergences in the VEVs of the local physical observables are well known in quantum field theory 
with boundaries and result from the idealization of the boundaries as perfectly smooth surfaces 
which are perfect reflectors at all frequencies. These divergences are investigated in detail for 
various types of fields and general shape of smooth boundary [20]l34j. Near the smooth boundary 
the leading divergence in the field squared varies as (D — l)th power of the distance from the 
boundary. It seems plausible that such effects as surface roughness, or the microstructure of the 
boundary on small scales (the atomic nature of matter for the case of the electromagnetic field 
[35] ) can introduce a physical cutoff needed to produce finite values of surface quantities. 

Now we turn to the VEVs of the energy-momentum tensor. By making use of formula ([3]), 
for the non-zero components one obtains (no summation over i) 

(Ti)M r(gj*) f 1M gM)rg) ^ cos^ 

{ 1 ^ " 2 ^%^+i U U sin^ 1 # ' { 2 ;ren 2V^r" s in D ^ ' 

(24) 

where i = 0, 1, . . . , D, and we use the following notations 



f® = 1 + (4£ - 1) [(D - l)S dl sin 2 + D5 J2 \ , % = 0, 3, . . . , D, (25) 

$ = 4 0) - ~ &) sin2 ^ > fj? = D [ i sin2 <i>? & ~ W - • ( 26 ) 

In the case 4>o = tt/2 and for minimally and confer mally coupled scalar fields, it can be checked 
that from formulae (|24p . after the transformation from cylindrical coordinates to the cartesian 



,2 Jj) 
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ones, as a special case we obtain the result derived in For a conformally coupled scalar field 



(0 

2/ 



and from (1241) one finds 



(T, 



k\(w) 
i I ren 



2£>+2 7r nr i r D+i 



L>-(L>-l)sin 2 (7rV^) 
L>sin D+1 (vr//g) 



diag(l, !,-£>,!, 



(27) 



In this case the vacuum energy-momentum tensor does not depend on the angular coordinate. 
For a non-conformally coupled field the VEVs (|24p diverge on the boundaries 4> = <p m and for 
points away from the edge r = 0, these divergences are the same as those for the geometry of a 
single plate. 

In the most important case D = 3, for the components of the renormalized energy-momentum 
tensor we find 



(T ° 



0\O) 

ren 



\ 1 /ren 



\ 1 2 /ren 



\ 1 2 /ren 



1 



3 /ren 32^ 2 r 4 



1-Q 4 
45 



+ 12 



1 



sin 2 (q<p) 



1 



32vr 2 r 4 ^ 45 

(C-Cc)g 2 



sin (q<f>) 
--U l -9 2 )(C-Cc) 



2 2 2 
3* + 3 



+ 12 



sin (q(p) 



1 



sin ((70) 



3(£-&) 9 3 cos (#) 



8vr 2 r 3 



sm 



8vr 2 r 4 



1 



60 



+ (e - e c ) i - ^ + 



3g 2 



sin (g^) 



(28) 



(29) 
(30) 
(31) 



Though we have derived these formulae for integer values of the parameter q, by the analytic 
continuation they are valid for non-integer values of this parameter as well. For a conformally 
coupled scalar field we obtain the result previously derived in the literature |19l I20j . The 
corresponding vacuum forces acting on the wedge sides are determined by the effective pressure 
— (T^ien- These forces are attractive for the wedge with q > 1 and are repulsive for q < 1. 



4 Field squared and the energy-momentum tensor 

We now turn to the geometry of a wedge with additional cylindrical boundary of radius a. 
Taking the coincidence limit x' — > x in formula (llip for the Wightman function and integrating 
over k, the VEV of the field squared is presented as the sum of two terms: 

(0|v? 2 |0) = (O^ 2 ^) + (v? 2 ) cyl , (32) 
where the part induced by the cylindrical boundary is given by the formula 

2 2^7^^ . 2 f°° 2 2 ^ K qn (az) t2 

\<P )cyi = - , D _ n ~ >sm (gn0) / dzz(z -m) I qn \ rz )- (33) 

Note that this part vanishes at the wedge sides <j> = 4> m , ^ r < a. Near the edge r = the 
main contribution into (( / 9 2 ) cy i comes from the term n = 1 and (^ 2 )cyl behaves like r 2q . The part 
{<f 2 ) cy i diverges on the cylindrical surface r = a. Near this surface the main contribution into 
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comes from large values n and for | (ft — <ft m | 3> 1 — r/a the leading behavior is the same as 
that for a cylindrical surface of radius a. 

Similarly, the VEV of the energy-momentum tensor for the situation when the cylindrical 
boundary is present is written in the form 

(0|T tt |0) = (0 w \T lk \0 w ) + (T ik ) cyh (34) 

where {Tik) cy \ is induced by the cylindrical boundary. This term is obtained from the corre- 
sponding part in the Wightman function, ((p(x)(p(x')) cy i, by using formula ([3]). For points away 
from the cylindrical surface this limit gives a finite result. For the corresponding components of 
the energy-momentum tensor one obtains (no summation over i) 

(47r) -V - /•«>,, ,^K qn {az) 



1 ( — j <fto n=l Jm 



Iqn(az) 

x {«£itV(^)] " cos(2gn^)} , (35) 

, mh 2(47t) _:2 2 Jl ^ . f°° , 9 9n d-3 K an (az) 



xl qn (rz) 



2 ll qn {rz) + {l-Ai)l' qn (rz) 



(36) 



with the notations 



a^[g(y)] = (^-l)[g' 2 (y) + (l±l 2 /y 2 )g 2 (y)]+2g 2 (y) 1 mV/y2 



D-1 7 

= 5 ' 2 (y) + (4e/yMyy(y)-5 2 (y){i±[i-4e(iTi)]/ 2 /y 2 }, (37) 

ag } b(y)] = (4£ - 1) [ 5 ' 2 (y) + g 2 (y)] - {At/y)g{y)g> \y) + g 2 {y) (4£ ± 1) /7y 2 , 

for a given function i = 0, 3, ...,D. In accordance with the problem symmetry, the 

expressions for the diagonal components are invariant under the replacement (ft — > 0o — (ft, and 
the off-diagonal component (7^ ) cy i changes the sign under this replacement. Note that the latter 
vanishes on the wedge sides (ft = (ft m , ^ r < a and for (ft = (fto/2. On the wedge sides for the 
diagonal components of the energy-momentum tensor we obtain (no summation over i) 

^W^ = ^S^En 2 rdzz(z 2 -rn 2 )^^^I 2 qn (rz), (38) 

1 ^ — 2 — / ^ ^0 n=l Jm 1 qn{UZ) 

where Ai = 4£ — 1, i = 0, 1, 3, . . . , D, A2 = 1. In particular, the additional vacuum effective 
pressure in the direction perpendicular to the wedge sides, p a = — (^2 )cyi.0=</> m j does not depend 
on the curvature coupling parameter and is negative for all values < r < a. This means that 
the vacuum forces acting on the wedge sides due to the presence of the cylindrical boundary are 
attractive. The corresponding vacuum stresses in the directions parallel to the wedge sides are 
isotropic and the energy density is negative for both minimally and conformally coupled scalars. 

For < r < a the cylindrical parts (|35p and (136j) are finite for all values ^ (ft ^ (fto, 
including the wedge sides. The divergences on these sides are included in the first term on the 
right-hand side of (|34p corresponding to the case without cylindrical boundary. Near the edge 
r = the main contribution into the boundary parts comes from the summand with n = 1 and 
one has (T-) cy \ oc r 2g ~ 2 , (T 2 1 ) cy i oc r 2q ~ x . The boundary part (2$ ) j diverges on the cylindrical 
surface r = a. Expanding over a — r, on the wedge sides for the diagonal components one finds 

A T (^ ±1 ) 

(g?W=*»* — »V. 2 > r-o, (39) 

2(47r)^~(a - r) D+l 
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where the coefficients are defined in the paragraph after formula (|38|) . It can be seen that 
for the off-diagonal component to the leading order one has (T 2 1 ) a oc (a — r)~ D . For angles 
< (j) < (fto, and for \<j>— <j> m \ 3> 1 — r/a, the leading divergence coincides with the corresponding 
one for a cylindrical surface of the radius a. 

Taking the coincidence limit of the arguments, from formula (|17p we obtain the VEV of the 
field squared in the region r > a: 

2 3-D °° f oo D _ 3 j / \ 

(^)cyl = - i (D1 , A E sin2 (^) / dzz (** ~ m ")^~ ¥f\ K Urz\ (40) 

As for the interior region, the VEV (|40p diverges on the cylindrical surface. For large distances 
from the cylindrical surface, r 3> a, and for a massless field the main contribution comes from 
the n = 1 term and to the leading order one finds (y 2 ) C yi ^ (a/r) D ~ 1+2q . For a massive field 
and for mr ^> 1 the part (^ 2 ) cy i is exponentially suppressed. 

For the part in the vacuum energy-momentum tensor induced by the cylindrical surface in 
the region r > a, from (J3]), (|17p . (|40p one has the following formulae 



x {«S^[^n(rz)] - a^i^rz)] cos(2gn0)} , (41) 
(Tl )cyl = f>sin(2 g n0) /°° dz z\z* - rn^ ^± 

1 (— 5— J 00 n=1 Jm i^qn\aZ) 



xK qn (rz) 



^K qn (rz) + (l-^)K' qn (rz) 



(42) 



with the functions ^^[^(y)] defined by (|37|). In the way similar to that used above for the VEV 
of the field square, it can be seen that at large distances from the cylindrical surface, r S> a, 
the main contribution comes from the term with n = 1 and for a massless field the components 
of the induced energy-momentum tensor behave as (T?) cy i oc (a/r) D+1+2q , (TrJ ) cy i oc (a/r) D+2q . 
As for the interior region, the vacuum forces acting on the wedge sides due to the presence of 
the cylindrical shell are attractive and the corresponding energy density is negative for both 
minimally and conformally coupled scalar s. 

In the limit (fio —* 0, r, a — > oo, assuming that a — r and a0o = b are fixed, from the results 
given above we obtain the vacuum densities for the geometry of two parallel plates separated by 
a distance b, perpendicularly intersected by the third plate. The vacuum expectation values of 
the energy-momentum tensor for this geometry of boundaries are investigated in [36] for special 
cases of minimally and conformally coupled massless scalar fields. 



5 VEVs for the electromagnetic fields 
5.1 Interior region 

In this section we consider a wedge with a coaxial cylindrical boundary assuming that all bound- 
aries are perfectly conducting. For this geometry there are two different types of the eigenfunc- 
tions corresponding to the transverse magnetic (TM, A = 0) and transverse electric (TE, A = 1) 
waves. In the Coulomb gauge, the vector potentials for these modes are given by the formulae 

a -a / ( l l iuJ ) (7 2e 3 + ikVt) J qn hr)sm(qn<p)ex.p [i(kz - ut)) , A = . . 

a la \ -e 3 xV t {J qn br)cos(qncj>)exp[i(kz-u;t)}}, A = l ' ^ 
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where 63 is the unit vector along the axis of the wedge, Vt is the part of the nabla operator 
transverse to this axis, and u> 2 = j 2 + k 2 . In Eq. (|43p . n = 1, 2, . . . for A = and n = 0, 1, 2, . . . 
for A = 1. From the normalization condition one finds 



0, 



1/2, 
1, 



n = 
n / 



jj^) + (l-z,2/x 2 )J 2 (*) 



(44) 
Eigenfunc- 



where we have introduced the notation T v (x) = x 

tions ()43p satisfy the standard boundary conditions on the wedge sides. From the boundary 
conditions on the cylindrical shell it follows that the eigenvalues for 7 are roots of the equation 



where jj) \x) 



jW( 7 a) = 0, A = 0,1, (45) 
= J' u {x). We will denote the corresponding eigenmodes by 



= J u (x) and Jv {x) 
7« = A^,J = 1,2,.... 

First we consider the VEVs of the squares of the electric and magnetic fields inside the shell 
Substituting the eigenfunctions (|43p into the corresponding mode-sum formula, we find 



<0|F 2 |0> 



m=0 



E *EE 



00 \(^)3m 



qui 



A=0,1 n=l 



^*$tt),J qn (\^r/a)], (46) 



where F = E,B with ijex = A, t]b\ = 1 — A, and the prime in the summation over n means that 
the term n = should be halved. In formula (1461) we have introduced the notations 



<7< o >[*(0),/(s)] 
g^{<j>)J{x)\ 



and 



(A;V/x 2 ) [$ 2 (0)/ /2 (x) + <D' 2 (0)/ 2 (x)/x 2 ] + ^ 2 (0)/ 2 (x) 
(1 + k 2 r 2 /x 2 ) [$ 2 (<P)f' 2 (x) + $ /2 (0)/ 2 (x)/x 2 ] , 

$W(<f>) = sin(z^), $W(<^) = cos(i/0). 



(47) 



(48) 



The expressions (|46h corresponding to the electric and magnetic fields are divergent. They may 
be regularized introducing a cutoff function ip^(uj) with the cutting parameter [i which makes 
the divergent expressions finite and satisfies the condition i/Jn(u)) — > 1 for fi — > 0. After the 
renormalization the cutoff function is removed by taking the limit fi — » 0. 

In order to further simplify the VEVs, we apply to the series over n the summation formula 
(jlOP for the modes with A = and the similar formula from [33] for the modes with A = 1. As 
it can be seen, for points away from the shell the contribution to the VEVs coming from the 
second integral terms on the right-hand sides of these formulae are finite in the limit fi — > 
and, hence, the cutoff function in these terms can be safely removed. As a result the VEVs are 
written in the form 

(0|F 2 |0) = (0 W \F 2 \0 W ) + (F 2 " 



'cyl 



(49) 



where 
(0 W 



F 2 |0„ 



E 

n=0 



dk 



dr.- 



vV + k 2 



1 + 



2^ 
T 2 " 



+JgnM ~ (- 1 )' ?F1 cos(2gn0) 



1 + 



q 2 n 2 

sy 2 



JanM 



(50) 



and 



= -T! E r dxx^^^G^\^{<P)J qn {xr)]. (51) 
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In formula (|51|) we have introduced the notations 

G<°> /(a:)] = $\<t>)f'\x) + f\x)/x 2 + 2<f> 2 (<t>)f(x), 

aWf*^),/^)] = -$ 2 (^)/' 2 (x)-ci>' 2 ^)/ 2 (x)/ a; 2 . (52) 

The second term on the right-hand side of Eq. (|49[) vanishes in the limit a — > oo and the first 
one does not depend on a. Thus, we can conclude that the term (0 W \F 2 \0 W ) corresponds to the 
part in the VEVs when the cylindrical shell is absent. 

First, let us concentrate on the part corresponding to the wedge without a cylindrical shell. 
In (|50j) the part which does not depend on the angular coordinate eft is the same as in the 
corresponding problem of the cosmic string geometry with the angle deficit 2ir — 4>q (see |32|). 
which we will denote by (0 S |-F 2 |0 S ). For this part we have 



<0 S 



|F 2 |0 S 



(OmI^IOm) 



(g 2 -l)(g 2 + ll) 
180vrr 4 



(53) 



where (0m|-F 2 |0m) is the VEV in the Minkowski spacetime without boundaries and in the last 
expression we have removed the cutoff. To evaluate the part in (|50p which depends on <p, we 
firstly consider the case when the parameter q is an integer. In this case, the summation over 
n can be done explicitly and the integrals are evaluated by introducing polar coordinates in the 
(k, 7)-plane. As a result, for the renormalised VEVs of the field squared in the geometry of a 
wedge without a cylindrical boundary we find 



(F 



2\H 



(g 2 -l)(g 2 + ll) 
180vrr 4 



27rr 4 sin 2 (g0) 



1 2 + 



3q 2 



2 sur (<?(/>) 



(54) 



with t)e\ = 1 and tjbi = 0. Though we have derived this formula for integer values of the 
parameter q, by the analytic continuation it is valid for non-integer values of this parameter as 
well. The expression on the right of formula (|54f) is invariant under the replacement eft — > (j)Q — <f> 
and, as we could expect, the VEVs are symmetric with respect to the half-plane (f> = (f>o/2. 
Formula (1541) for F = E was derived in Ref. [22] within the framework of Schwinger's source 
theory. 

Now, we turn to the investigation of the parts in the VEVs of the field squared induced by 
the cylindrical boundary and given by formula (|5ip . These parts are symmetric with respect to 
the half-plane 4> = <Ao/2- The expression in the right-hand side of (|5ip is finite for < r < a 
including the points on the wedge sides, and diverges on the shell. To find the leading term 
in the corresponding asymptotic expansion, we note that near the shell the main contribution 
comes from large values of n. By using the uniform asymptotic expansions of the modified 
Bessel functions for large values of the order, up to the leading order, for the points a — r <C 
a| sin0|, a\ sin(<^o — 4>)\ we find {F 2 ) cy \ ~ — 3(— l) VF1 /[An(a — r) 4 ]. These surface divergences 
originate in the unphysical nature of perfect conductor boundary conditions. In reality the 
expectation values will attain a limiting value on the conductor surface, which will depend on 
the molecular details of the conductor. From the formulae given above it follows that the main 
contribution to (F 2 ) cy i are due to the frequencies u < (a — r) _1 . Hence, we expect that formula 
(f5"T|) is valid for real conductors up to distances r for which (o — r) _1 -C u>o, with luq being the 
characteristic frequency, such that for oj > ujq the conditions for perfect conductivity fail. 

Near the edge r = 0, assuming that r/a ^ 1, the asymptotic behavior of the part induced 
in the VEVs of the field squared by the cylindrical shell depends on the parameter q. For 
q > 1 + rjpi, the dominant contribution comes from the lowest mode n = and to the leading 
order one has (F 2 ) cy \ oc r 2r?F1 . In this case the quantity (B 2 ) cy \ takes a finite limiting value 
on the edge r = 0, whereas (E 2 ) cy \ vanishes as r 2 . For q < 1 + rjpi the main contribution 
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comes from the mode with n = 1 and the shell-induced parts diverge on the edge r = with 
(F 2 ) cy \ cc r 2 ^ -1 ). In accordance with (|54p . near the edge r = the total VEV is dominated by 
the part coming from the wedge without the cylindrical shell. 



5.2 Exterior region 



In the exterior region (region II in figure [T|) , the corresponding eigenfunctions for the vector 
potential are obtained from formulae (|43p by the replacement 



Jgndr) -> s$( 7 a,7r) = J qn (ir)Yg\ ia ) - Y qn ( ir )J$ ( 7 a) 



(55) 



where, as before, A = 0, 1 correspond to the waves of the electric and magnetic types, respectively. 
The eigenvalues for 7 are continuous and 



(3- 2 = (87T/q)S nl u \&{ia)+Y$ 2 { ia ) 



(56) 



Substituting the eigenfunctions into the corresponding mode-sum formula, for the VEV of the 
field squared one finds 



<0|F 2 |0> 



t^J-oo Jo 7 A t^\/FTy Jg )2 (ja)+Yff 2 ^a) 



n=0 



(57) 



where the functions g^ FX ^ [$(</>), f(x)] are defined by relations (|4T|) with /(x) = 5^(70, x). To 
extract from this VEV the part induced by the cylindrical shell, we subtract from the right- 
hand side the corresponding expression for the wedge without the cylindrical boundary. As a 
result, the VEV of the field squared is written in the form (|49j) . where the part induced by the 
cylindrical shell is given by the formula 



{F 2 



>cyl 



2? W 

7T ^ 



E 

A=0,1 



dxx zAn{xa 



K$ (xa) 



q(vfx) [$C\) (0) f K qn (xr)]. 



(58) 



In this formula the functions G^ r,FX ^ [$(0), f(x)] are defined by expressions (|52|) . Comparing this 
result with formula (|5ip . we see that the expressions for the shell-induced parts in the interior 
and exterior regions are related by the interchange I qn T± K qn . 

The VEV (|58p diverges on the cylindrical shell with the leading term being the same as 
that for the interior region. At large distances from the cylindrical shell we introduce a new 
integration variable y = xr and expand the integrand over a/r. For q > 1 the main contribution 
comes from the lowest mode n = and up to the leading order we have 



(E 2 



>cyl 



4q (a/rY 



(B 2 



57rr 4 ' x cyl 157rr 4 

For q < 1 the dominant contribution into the VEVs at large distances is due to the mode n = 1 
with the leading term 



28q (a/rf 



(59) 



irr H 



cos(2q(p) 
2q + 3 



9 + 1 
2q + l 



(60) 



For the case q = 1 the contributions of the modes n = and n = 1 are of the same order 
and the corresponding leading terms are obtained by summing these contributions. The latter 
are given by the right-hand sides of formulae (f59|) and ([601) . As we see, at large distances the 
part induced by the cylindrical shell is suppressed with respect to the part corresponding to the 
wedge without the shell by the factor (a/r) 213 with j3 = min(l,g). 
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6 Energy-momentum tensor for the electromagnetic field 

Now let us consider the VEV of the energy-momentum tensor in the region inside the cylindrical 
shell. Substituting the eigenfunctions (|43|) into the corresponding mode-sum formula, for the 
non-zero components we obtain (no summation over i) 

00 r+00 00 \( A ) 3 7^ (\( x ^\ 

a n=0 J -00 A=(J)1 j=1 ^ A W2 + k2a2 



00 



n=0 " * A=0,1 

A n,j I qn{^n,j) j2 / \ (A) 



X 2J / (X)2 W r/a), (62) 



"J 

where i = 0, 1, 2, 3, and we have introduced the notations 

/(x)] = (-l)*( 2 £; 2 / 7 2 + l) [$> 2 (<t>)f' 2 (x) + &\<t>)f 2 {x)/y 2 ] + d> 2 («£)/ 2 (x), 

= (-iy$ 2 (^)/ ,2 (x) - [$ 2 (<£) + (-iy$' 2 (^)/x 2 ] / 2 (x), (63) 

with j = 0, 3 and / = 1,2. As in the case of the field squared, in formulae (|6ip and (|62|) we 
introduce a cutoff function and apply formula (|10|) for the summation over n. This enables us 
to present the vacuum energy-momentum tensor in the form of the sum 

(0\T?\0) = (0 w \T t k \0 w ) + (T?)^, (64) 

where (0 M |T^|0 to ) is the part corresponding to the geometry of a wedge without a cylindrical 
boundary and {T^) cy \ is induced by the cylindrical shell. The latter may be written in the form 
(no summation over i) 

00 r°° \ 

(^)cyi = AE'E / ^x 3 %^F«[cI>W(0),V(xr)], (65) 

(xa) 



n=0 A=0,1 ' 



(^)cyi = ^^E'— ( 2 ^)E(-D A / ^%rr^ : 

47F n=0 A=0,1 J ° Iqn{xa) 



xr), (66) 



with the notations 

F® [*($,/(!/)] = <*> 2 (</>)/ 2 (y), ^ = 0,3, 

F® [$(«£), /(*/)] = -(-l)^ 2 (0)/«(y)-[^)-(-l)^(^)/i/ 2 ]/ 2 (y), t = l,2. (67) 

The diagonal components are symmetric with respect to the half-plane <fi = <fio/2, whereas the 
off-diagonal component is an odd function under the replacement cj> — > <j)Q — <fi. As it can be 
easily checked, the tensor (Zl ) C yl is traceless and satisfies the covariant continuity equation. The 
off-diagonal component (T 2 1 ) cy i vanishes at the wedge sides and for these points the VEV of the 
energy-momentum tensor is diagonal. The vacuum energy density induced by the cylindrical 
shell in the interior region is always negative. 

The renormalized VEV of the energy density for the geometry without the cylindrical shell 
is obtained by using the corresponding formulae for the field squared. Other components are 
found from the tracelessness condition and the continuity equation and one has [2], HHJ [20] 

(T*)l^ = - (g2 ~72oX^ 11) diag(1 ' _3 ' 1} ' (68) 
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Formula (|68|) coincides with the corresponding result for the geometry of the cosmic string with 
the angle deficit 2tt — (fro and in the corresponding formula q = 2tt/4>q. 

The normal force acting on the wedge sides is determined by the component (T| ) rcn of the 
vacuum energy- momentum tensor evaluated for = and <fr = (fro. On the base of formula (|64p 
for the corresponding effective pressure one has 

P2 = -(72 2 ) ren U=0,(/>o = P2w + £>2cyl, (69) 

where p2 W = — (T^ren is the normal force acting per unit surface of the wedge for the case 
without a cylindrical boundary and the additional term 

P 2cy i = -(T 2 2 ) cyl U =0 ,, =-^E E / dxx 3 ^^F^[I qn (xr)}, (70) 

n=0 A=0,1 Iqn \ Xa ) 

with the notations 

F^[f(y)] = v 2 f(y)/y 2 , F^[f( y )] = -f' 2 ( y ) - f\ y ), (71) 

is induced by the cylindrical shell. Note that the normal force on the wedge sides is the sum of the 
corresponding forces for Dirichlet and Neumann scalars corresponding to the terms with A = 
and A = 1 respectively. The finiteness of the normal stress on the wedge sides is a consequence 
of the fact that for a single perfectly conducting plane boundary this stress vanishes. Note 
that this result can be directly obtained from the symmetry of the corresponding problem with 
combination of the continuity equation for the energy-momentum tensor. It also survives for 
more realistic models of the plane boundary (see, for instance, |37| ) though the corresponding 
energy density and parallel stresses no longer vanish. So we expect that the obtained formula for 
the normal force acting on the wedge sides will correctly approximate the corresponding results 
of more realistic models in the perfectly conducting limit. The corresponding vacuum forces are 
attractive for q > 1 and repulsive for q < 1. In particular, the equilibrium position corresponding 
to the geometry of a single plate (q = 1) is unstable. As regards to the part induced by the 
cylindrical shell, from (f70l) it follows that P2cyi < and, hence, the corresponding forces are 
always attractive. In figure [2] we have plotted the vacuum pressure on the wedge sides induced 
by the cylindrical boundary versus r/a for Dirichlet scalar (left panel) and for the electromagnetic 
field (right panel). The full (dashed) curves correspond to the wedge with (fro = n/2 ((fro = Stt/2). 

Now, let us discuss the behavior of the boundary-induced part in the VEV of the energy- 
momentum tensor in the asymptotic regions of the parameters. Near the cylindrical shell the 
main contribution comes from large values of n and for the points a — r <C a\ sin <j>\ , a\ sin(^o — <fr) \ 
the leading terms are the same as those for a cylindrical shell when the wedge is absent. For 
points near the edges (r = o, (fr = 0, 0q) the leading terms in the corresponding asymptotic 
expansions are the same as for the geometry of a wedge with the opening angle (fro = tt/2. 
Near the edge, r — * 0, for the components (no summation over i) (T-) cy \, i = 0, 3, the main 
contribution comes from the mode n = and we find (T?) cy i ~ — 0.0590g/a 4 , i = 0,3. For 
the components (no summation over i) (T?) cy i, i = 1,2, when q > 1 the main contribution 
again comes form n = term and one has (T?) cy i ~ — (Tg)^, i = 1,2. For q < 1 the main 
contribution into the components {T-) cy \, i = 1,2, comes from the term n = 1 and we have 
(no summation over i) (T-) cy \ cc r 2 ^ -1 ), i = 1,2. In this case the radial and azimuthal stresses 
induced by the cylindrical shell diverge on the edge r = 0. For the off-diagonal component the 
main contribution comes from the n = 1 mode and (T^cyi behaves like r 2q ~ l . 

Now we turn to the VEVs of the energy-momentum tensor in the exterior region. Subtracting 
from these VEVs the corresponding expression for the wedge without the cylindrical boundary, 
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Figure 2: The effective azimuthal pressure induced by the cylindrical shell on the wedge sides, 
° 4 P2cyh as a function of rja for Dirichlet scalar (left panel) and for the electromagnetic field 
(right panel). The full (dashed) curves correspond to q = 2 (g = 2/3). 

analogously to the case of the field square, it can be seen that the VEVs are presented in the 
form (I64p . with the parts induced by the cylindrical shell given by the formulae (no summation 
over i) 

(^)cyi = AE'E / dxx^^^F^^K^xr)), (72) 
Z7T „=o AM),!" 70 K K qn >(xa) 

2 £i 00 roo t-(A) / \ 

(Tl) cyl = ^E'nsin^) £(-!)*/ dxx ~S^ K U^)- (73) 
47r ar ™-n a=o,i Jo K K qn '{xa) 



m=0 



Here the functions [$(0), f(y)] are defined by expressions (|671) . It can be seen that the 
vacuum energy density induced by the cylindrical shell in the exterior region is positive. 

In the way similar to that for the interior region, the force acting on the wedge sides is 
presented in the form of the sum (|69p . where for the part due to the presence of the cylindrical 
shell we have 

p 2cy i = -<T 2 2 > cyl U =0 ,, = -4E E / dxx"^g^Fg)[K qn {xr)]. (74) 

n=o a=o,i Jo K y qn '{xa) 

In this formula, the function F^ [/(?/)] is defined by relations (|7ip and the corresponding forces 
are always attractive. 

The leading divergence in the boundary induced part (|72|) on the cylindrical surface is given 
by the same formulae as for the interior region. For large distances from the shell and for q > 1 
the main contribution into the VEVs of the diagonal components comes from the n = 0, A = 1 
term and one has (no summation over i) 

/rri\ _ («/ r ) 2 i r (k\ 

\Ti ) cy i ~ 157r2r4 , c = C 3 = 2, Cl = 1, c 2 = -5. (75) 

In the case q < 1 the main contribution into the VEVs of the diagonal components at large 
distances from the cylindrical shell comes from the n = 1 mode. The leading terms in the 
corresponding asymptotic expansions are given by the formulae 

/t*\ ~ 2/ , n , x cos(2g0) (a\^i 

{Ti )cyl « -q (q + l)ci(g) 2 4 ^- j , (76) 
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with the notations 



co (<?) = c 3 (q) 



1 



ci(<?) 



2q 2 + q + 1 



C2(<?) 



q+1 
2q+l 



(77) 



2g + 3 



(2g + l)(2 (? + 3) 



In the case q = 1 the asymptotic terms are determined by the sum of the contributions coming 
from the modes n = and n = 1. The latter are given by formulae (|75p . (|76j) . For the off- 
diagonal component, for all values q the main contribution at large distances comes from the 
n = 1 mode and (T 2 1 ) cyl oc (a/r) 2q r~ 3 . 

7 Conclusion 

We have investigated the polarization of the scalar and electromagnetic vacua by a wedge with 
coaxial cylindrical boundary, assuming Dirichlet boundary conditions in the case of a scalar field 
and perfectly conducting boundary conditions for the electromagnetic field. The application of 
the Abel-Plana-type summation formula for the series over the zeros of the Bessel function and its 
derivative allowed us to extract from the VEVs the parts due to the wedge without a cylindrical 
boundary and to present the additional parts induced by this boundary in terms of exponentially 
convergent integrals. The vacuum densities for the geometry of a wedge without a cylindrical 
boundary are considered in section [3l We have derived formulae for the renormalized VEVs 
of the field squared and the energy- momentum tensor, formulae (I22p . (|24p . For a conformally 
coupled scalar the energy-momentum tensor is diagonal and does not depend on the angular 
variable (ft. The corresponding vacuum forces acting on the wedge sides are attractive for (fto < tt 
and are repulsive for (fto > tt. 

For a scalar field the parts in the Wightman function induced by the cylindrical boundary 
are given by formulae f)13[) and (|17f) for the interior and exterior regions respectively. The 
corresponding VEVs for the field squared and the energy-momentum tensor are investigated in 
section [H The field squared is given by formula (133 [) and vanishes on the wedge sides (ft — (ft m 
for all points away from the cylindrical surface. The energy-momentum tensor induced by 
the cylindrical surface is non-diagonal and the corresponding components are determined by 
formulae ([35]) . ([36]) . The off-diagonal component vanishes on the wedge sides. The additional 
vacuum forces acting on the wedge sides due to the presence of the cylindrical surface are 
determined by the ^-component of the corresponding stress and are attractive for all values 
(fto- On the wedge sides the corresponding vacuum stresses in the directions parallel to the 
wedge sides are isotropic and the energy density is negative for both minimally and conformally 
coupled scalars. The formulae in the exterior region differ from the corresponding formulae for 
the interior region by the interchange I qn {z) ^ K qn (z). For large distances from the cylindrical 
surface, r 3> a, the VEVs behave as (a/r) D ~ 1+2q for the field squared and as (a/r) D+l+2q for 
the diagonal components of the energy-momentum tensor. 

In the second part of the paper we have evaluated the VEVs of the field squared and the 
energy-momentum tensor for the electromagnetic field. For the wedge without the cylindrical 
shell the VEVs of the field squared are presented in the form (|54p . The first term on the right 
of this formula corresponds to the VEVs for the geometry of a cosmic string with the angle 
deficit 2-7T — (J)q. The parts induced by the cylindrical shell are presented in the form (|5ip for 
the interior region and in the form (j58h for the exterior region. We have discussed these general 
formulae in various asymptotic regions of the parameters. In section [6] we consider the VEV of 
the energy-momentum tensor. For the geometry of a wedge without the cylindrical boundary 
the vacuum energy-momentum tensor does not depend on the angle (ft and is the same as in 
the geometry of the cosmic string. The corresponding vacuum forces acting on the wedge sides 
are attractive for (fto < tt and repulsive for (fto > tt. In particular, the equilibrium position 
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corresponding to the geometry of a single plate is unstable. For the region inside the shell the 
part in the VEV of the energy-momentum tensor induced by the presence of the cylindrical 
shell is non-diagonal and the corresponding components are given by formulae (|65|) . (|66|) for the 
interior region and by (|72p . (|73p for the exterior region. The vacuum energy density induced 
by the cylindrical shell is negative in the interior region and is positive in the exterior region. 
For a wedge with cpo < it the part in the vacuum energy-momentum tensor induced by the shell 
is finite on the edge r = 0. For 0o > vr the shell-induced parts in the energy density and the 
axial stress remain finite, whereas the radial and azimuthal stresses diverge as r 2 ( 7r /</ , o-i) i ^he 
corresponding off-diagonal component behaves like r 27T ^°~ 1 for all values 4>q. For points near 
the edges (r = a, eft = 0, 0o)> the leading terms in the corresponding asymptotic expansions are 
the same as for the geometry of a wedge with the opening angle <j)Q = tt/2. The presence of the 
shell leads to additional forces acting on the wedge sides. The corresponding effective azimuthal 
pressures are given by formulae (|70p . (|74p and these forces are always attractive. 
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